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[i]  Modulated  HE  heating  of  the  auroral  electrojet  is  used 
to  inject  ELF  signals  into  the  magnetosphere  that  are 
observed  on  the  low  altitude  DEMETER  spacecraft.  The 
HE  heater  is  a  component  of  the  High-Frequency  Active 
Auroral  Research  Program  (HAARP)  facility  in  Gakona, 
Alaska,  (located  at  L  ~  4.9).  Simultaneous  observations  of 
all  six  components  of  the  ELF  electromagnetic  fields  on  the 
DEMETER  spacecraft  are  used  to  estimate  the  total  ELF 
power  radiated  downward  into  the  Earth-ionosphere 
waveguide.  ELF  signals  generated  by  the  HAARP  heater 
are  also  simultaneously  observed  at  a  nearby  ground-based 
site,  allowing  a  comparison  of  the  ELF  power  in  the  Earth- 
ionosphere  waveguide  versus  that  detected  on  DEMETER. 
The  estimated  values  of  power  onboard  DEMETER  at 
different  frequencies  range  from  0.32W  to  4W,  while  the 
values  of  power  estimated  from  a  ground  receiver  at  a 
distance  of  36  km  from  HAARP  range  from  2.7  IW  to 
4.22W.  Citation:  Platino,  M.,  U.  S.  Inan,  T.  F.  Bell,  M.  Parrot, 
and  E.  J.  Kennedy  (2006),  DEMETER  observations  of  ELF  waves 
injected  with  the  HAARP  HE  transmitter,  Geophys.  Res.  Lett,  33, 
L 1 6 1 0 1 ,  doi:  1 0. 1 029/2006GL026462. 

1.  Introduction 

[2]  Generation  of  ELF/VLF  waves  by  modulating  the 
D  region  conductivity  in  the  vicinity  of  the  auroral  electrojet 
using  a  powerful  High  Frequency  (HE)  heater  is  by  a  now  well 
established  technique  [Stubbe  and  Kopka,  1977;  Rietveld  et 
al,  1989;  Inan  et  al,  2004,  and  references  therein].  The 
mechanism  involves  the  periodic  modification  of  the  electron 
temperature  in  the  D  and  lower  E  regions  of  the  ionosphere, 
which  leads  to  periodic  changes  in  the  Pedersen  and  Hall 
conductivities,  resulting  in  periodically  changing  electric 
current  density  in  the  electrojet  and  thus  radiation.  Under 
such  circumstances,  the  heated  ionospheric  area  with  its 
immediate  surroundings  acts  as  a  huge  ‘polar-electrojet ’ 
antenna,  radiating  primarily  at  a  frequency  corresponding  to 
the  modulation  frequency  of  the  HE  carrier. 

[3]  Spacecraft  observations  of  ELFA^LF  waves  generated 
by  different  HE  heating  facilities  have  been  carried  out  on 
both  low  [Lefeuvre  et  al,  1985;  James  et  al,  1984;  Inan  and 
Helliwell,  1985],  and  high-altitude  spacecraft  [James  et  al, 
1990;  Bell  et  al,  2004;  Platino  et  al,  2004].  Nevertheless, 
the  degree  to  which  modulated  HE  heating  can  be  consis- 
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tently  used  to  inject  ELF/VLF  waves  upward  into  the 
overlying  magnetosphere,  and  the  total  ELF  power  radiated 
by  such  a  source  is  not  understood  at  a  quantitative  level. 

[4]  In  this  paper,  we  present  recent  observations  on  the 
DEMETER  spacecraft  of  electromagnetic  ELF  waves  pro¬ 
duced  through  modulation  of  auroral  electrojet  currents,  and 
use  the  simultaneously  observed  multiple  field  components 
to  estimate  the  total  ELF  power  radiated  into  the  Earth- 
ionosphere  waveguide  by  the  electrojet  current  modified  by 
the  HAARP  HE  heater.  The  observations  reported  herein 
were  carried  out  over  the  ELF  frequency  range  of  1 1 8  Hz  to 
1197  Hz  on  the  3 -axis  stabilized  DEMETER  satellite,  at  an 
altitude  of  ^730  km  in  an  orbit  of  98.3°  of  inclination. 
Ground  observations  of  the  HAARP  ELF  transmitted  sig¬ 
nals  are  also  provided  and  compared  to  the  DEMETER 
observations.  A  numerical  raytracing  model  was  used  to 
simulate  the  propagation  of  ELF  whistler  mode  waves  from 
the  ionospheric  generation  regions  up  to  the  spacecraft  and 
to  estimate  the  spatial  extent  of  the  region  in  the  ionosphere 
illuminated  by  ELF  waves. 


2.  HAARP-DEMETER  Campaign 

[5]  The  DEMETER  spacecraft  is  a  low  altitude  satellite, 
with  two  sets  of  antennas,  one  set  to  measure  the  three 
components  of  the  plasma  wave  electric  field  using  the 
electric  sensors  of  the  Electric  Field  Instrument  (ICE)  and 
one  set  to  measure  the  three  components  of  the  plasma  wave 
magnetic  field  using  the  magnetic  sensors  of  the  Instrument 
Magnetometer  Search-Coil  (IMSC),  operating  over  a  fre¬ 
quency  range  of  10  Hz  to  1.25  kHz  in  the  burst  mode 
[Parrot,  2002].  Attempts  to  observe  the  HAARP  ELF 
signals  on  DEMETER  were  carried  out  during  passes  close 
to  the  HAARP  facility,  on  spacecraft  trajectories  whose 
subsatellite  points  on  the  ground  were  at  a  distance  no  larger 
than  200  km  from  HAARP  at  the  closest  point.  The 
HAARP  HF  heater  is  located  near  Gakona,  Alaska,  at  a 
site  with  geographic  coordinates  of  62.4°  (North)  latitude, 
145.2°  (West)  longitude  and  geomagnetic  coordinates  of 
63.1°  (North)  latitude  and  92.4°  (West)  longitude.  At 
HAARP,  an  HF  phased  -  array  transmitter  is  used  to  heat 
small,  well-defined  volumes  of  the  ionosphere  with  high 
power  radio  signals  in  the  2.8  to  10  MHz  frequency 
range.  At  the  time  these  experiments  were  performed,  the 
HAARP  HF  transmitter  had  a  total  maximum  radiated 
power  capability  of  960  kW.  The  ELF/VLF  signal  format 
shown  in  Figure  1  was  impressed  upon  a  3.3  MHz  carrier 
through  sinusoidal  amplitude  modulation,  with  a  modula¬ 
tion  index  of  0.95.  The  HF  carrier  frequency  was  chosen  to 
provide  maximum  heating  in  the  D  region  of  the  ionosphere 
[James  et  al,  1984],  taking  into  account  the  allocated 
frequencies  at  which  HAARP  can  transmit  at  its  highest 
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Figure  1.  Schematic  spectrogram  (to  scale)  of  the 
transmission  pattern  used  by  the  HAARP  HE  facility  to 
heat  the  ionosphere.  The  ELF  waves  were  AM  modulated 
onto  a  carrier  frequency  of  3.3  MHz.  The  transmitted  HE 
wave  was  in  the  X  mode. 

power,  using  HE  waves  polarized  in  the  extraordinary 
(X)  mode  to  provide  the  maximum  ELF  radiated  power 
{Kapustin  et  ah,  1977;  Stubbe  et  al,  1981,  1982].  The  HE 
beam  was  oriented  vertically  in  order  to  have  a  transmitted 
HE  wave  with  a  k- vector  as  perpendicular  as  possible  to  the 
D  layer  of  the  ionosphere.  The  modulation  pattern  consists 
of  pulses  of  2  sec,  4  sec,  5  sec  and  12  sec,  duration  at 


frequencies  located  in  the  ELF  frequency  range.  These  ELF 
modulation  frequencies  are  below  the  range  of  frequencies 
which  were  found  in  the  past  to  be  more  efficiently 
generated  via  modulated  HE  heating  {Stubbe  et  al,  1981], 
but  all  lie  in  the  range  for  which  DEMETER  can  record 
in  the  burst  mode  the  three  components  of  both  electro¬ 
magnetic  fields  for  frequencies  below  1 .25  kHz.  The  actual 
numerical  values  were  chosen  to  avoid  possible  interference 
from  harmonics  of  60  Hz  generated  by  power  grids 
{Helliwell  et  al,  1975]. 

3.  Determination  of  the  Illuminated  Region 

[6]  The  raytracing  program  used  for  the  computations 
presented  here  is  described  by  Inan  and  Bell  [1977],  with  a 
geomagnetic  field  model  based  on  a  centered  dipole  with 
electron  gyrofrequency  of  880  kHz  at  the  ground  on  the 
magnetic  equator.  The  model  of  the  magneto  spheric  elec¬ 
tron  and  ion  densities  is  the  same  as  used  by  Platino  et  al 
[2004].  Nevertheless,  ray  paths  at  these  low  altitudes  are  not 
was  not  very  sensitive  to  the  magnetospheric  electron 
density  model. 

[v]  Figure  2  shows  a  schematic  of  the  DEMETER  pass 
over  HAARP.  The  propagation  path  of  the  wave  energy  is 
calculated  using  raytracing,  as  described  in  the  previous 
paragraph.  The  region  of  space  illuminated  by  HAARP  is 
determined  from  the  maximum  distance  at  which  the 
HAARP  signal  is  observed  on  the  spacecraft  (430  km  lateral 
distance  from  HAARP).  DEMETER  receivers  are  routinely 
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Figure  2.  Ray  tracing  simulation  of  the  wave  propagation  path  for  the  date  studied.  Shown  are  the  ray  paths  from  the 
simulation.  The  raypaths  were  injected  into  the  magnetosphere  at  very  low  wave  normal  angles  and  the  difference  in  paths 
for  the  different  frequencies  observed  was  not  noticeable  at  these  low  attitudes.  Included  in  the  plot  are  a  scheme  of  the 
satellite  path  during  February  10,  2005  at  around  06:52:30  UT,  and  different  views  of  the  illuminated  regions  from  the 
ground  in  the  bottom  left  panel  and  from  space  in  the  bottom  right  panel. 
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shut  off  at  invariant  latitudes  larger  than  65°,  so  that  no 
observations  exist  beyond  that  latitude.  We  assume  azi¬ 
muthal  symmetry  around  the  magnetic  field  line  intersecting 
HAARP  to  estimate  the  area  A  of  the  region  of  space 
illuminated  by  HAARP.  Using  backwards  ray  tracing  we  can 
then  determine  the  area  B  of  the  D-region  in  the  ionosphere 
illuminated  by  HAARP.  We  also  know  the  region  C  in  the 
ionosphere  heated  by  the  HE  transmitter,  using  the  known 
beamwidth  of  the  HAARP  array.  Since  ELF  waves  tend  to 
propagate  predominantly  along  the  Earth’s  magnetic  field 
lines,  upward  propagating  ELF  radiation  produced  in  and 
near  the  heated  spot  will  not  spread  significantly  in  the 
horizontal  direction  as  it  propagates  to  the  DEMETER 
altitude.  The  large  difference  between  regions  B  and  C,  is 
due  to  the  fact  that  waves  generated  by  the  modulated 
electrojet  are  produced  over  an  altitude  range  which  spans 
the  lower  ionospheric  reflection  height  (^85  km  at  night) 
and  not  only  propagate  upward,  but  also  downward.  The 
downward  propagating  waves  are  then  reflected  from  the 
ground  and  escape  into  the  magnetosphere  as  they  propa¬ 
gate  in  the  Earth-ionosphere  waveguide  to  larger  lateral 
distances  away  from  the  source.  This  effect  produces  an 
ELF  illumination  region  which  is  much  larger  than  the 
actual  HF-heated  ionospheric  spot.  Modulation  of  the  iono¬ 
spheric  conductivity  in  the  presence  of  the  auroral  electrojet 
electric  fields  can  also  lead  to  horizontal  currents  which 
flow  in  loops  outside  of  region  C,  producing  an  effective 
radiating  region  that  will  be  larger  than  the  actual  heated 
one.  To  take  into  account  these  looping  currents  requires  a 
model  of  the  entire  heating  process,  which  is  beyond  the 
scope  of  the  present  discussion. 
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Figure  3.  Detailed  ICE  and  IMSC  spectrograms  from  the 
DEMETER  satellite  pass  during  February  10,  2005.  (a)  An 
overview  of  the  detected  HAARP  signal  in  one  of  the 
electric  field  components.  The  transmitted  pattern  at 
frequencies  up  to  798  Hz  is  observed,  (b)  An  overview  of 
the  detected  HAARP  signal  in  one  of  the  magnetic  field 
components.  The  transmitted  pattern  at  frequencies  up  to 
798  Hz  is  observed. 


4.  Observations:  HAARP-Generated  ELF 
Signals  Observed  at  DEMETER  and  on 
the  Ground 

[s]  ELF/VLF  waves  generated  by  HAARP  -  modulated 
electrojet  currents  were  observed  on  one  fourth  of  the 
DEMETER  satellite  passes  over  HAARP  (3  times  of  12 
passes).  This  low  probability  may  be  due  to  the  variability 
of  the  intensity  of  the  electrojet  [Lefeuvre  et  al,  1985; 
James  et  al,  1984;  Inan  and  Helliwell,  1985],  as  well  as  the 
ambient  ELF  noise  levels  at  the  DEMETER  orbit.  However 
we  also  have  to  consider  the  fact  that  DEMETER  does  not 
record  data  above  the  invariant  latitude  of  65°,  which,  given 
the  fact  that  HAARP  is  located  at  62.4°  and  that  the  ELF 
raypaths  have  the  tendency  to  move  towards  higher  lat¬ 
itudes,  may  significantly  limit  the  chances  of  observing  the 
HAARP-generated  ELF  signal  to  those  occasions  when  the 
emitted  ELF  fields  are  relatively  intense. 

[9]  Figure  3  shows  an  overview  spectrogram  of  data 
recorded  using  the  ICE  and  IMSC  instruments  on 
DEMETER  during  February  10,  2005,  where  the 
HAARP-transmitted  signal  (Figure  I)  is  clearly  observed 
during  a  30  second  period.  Pulses  at  all  of  the  ELF 
modulation  frequencies  are  observed  except  the  one  at 
647  Hz  and  the  ones  above  800  Hz.  Between  600  Hz  and 
650  Hz,  a  strong  noise  band  is  seen  which  corresponds  to 
the  local  proton  gyrofrequency,  representing  the  upper  band 
limit  for  left  hand  mode  propagation.  The  647  Hz  signal 
from  HAARP  may  be  strongly  absorbed  (attenuated)  due  to 
its  proximity  to  the  upper  cutoff  since  the  647  Hz  signal  is 


observed  on  the  ground  at  the  same  levels  as  the  other  ELF 
pulses.  Pulses  at  frequencies  higher  than  800  Hz  are  not 
observed  because  the  DEMETER  receiver  stopped  record¬ 
ing  at  06:53:27  UT.  These  pulses  were  in  fact  also  produced 
and  launched  downward,  as  evidenced  by  the  fact  that  they 
are  observed  on  the  ground.  Similar  pulses  in  this  frequency 
range  were  actually  detected  on  DEMETER  during  other 
passes,  where  the  termination  of  the  data  recording  (at  the 
65 -degree  invariant  latitude  point)  was  at  a  later  time  in  the 
course  of  the  transmission  sequence.  It  can  be  seen  that 
many  pulses  exhibit  deep  fading,  especially  those  at  707  Hz 
and  798  Hz.  This  fading  is  possibly  due  to  multi -path 
effects  produced  by  ionospheric  irregularities  [Sonwalkar  et 
ai,  1984].  Natural  impulsive  emissions  are  observed  in  the 
electric  field  instrument  below  600  Hz.  A  resonance  at 
the  proton  gyrofrequency  is  observed  between  600  Hz  and 
650  Hz.  Other  artificial  emissions  are  observed  like  contin¬ 
uous  narrow  band  observed  around  180  Hz,  most  likely 
produced  by  power  line  harmonics.  The  spectral  lines 
observed  above  800  Hz  and  the  narrow  band  noise  observed 
below  50  Hz  in  the  magnetic  field  instrument  are  due  to 
instrumental  noise  onboard  the  satellite. 

[10]  Figure  4a  shows  simultaneous  data  recorded  on  the 
ground  near  Chistochina,  Alaska,  at  a  distance  of  ~36  km 
from  the  HAARP  HF  heater.  These  data  also  show  ELF 
signals  generated  by  HAARP.  Displayed  in  this  figure  is  the 
wave  magnetic  field  detected  with  one  of  the  two  loop 
antennas  on  the  ground.  This  type  of  data  provides  an 
indication  of  the  strength  of  the  ELFA^LF  signal  radiated 
downward  by  the  auroral-electrojet  antenna.  Since  the  extent 
of  the  generation  region  is  believed  to  encompass  the  upper 
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Figure  4.  (a)  Ground  observations  of  the  ELF  HAARP 
generated  signal  at  Chistochina,  on  February  10,  2005, 
starting  at  06:52:30  UT.  The  panel  shows  data  from  a 
vertical  magnetic  loop  antenna,  oriented  to  the  geographic 
East- West  line.  The  magnitude  scale  refers  to  magnetic  field 
intensity  in  pT  /  (b)  Scheme  used  for  the  calculation 

of  the  radiated  power  toward  the  ground.  Indicated  are  the 
electrojet  currents,  moving  outside  the  page,  and  the 
location  of  the  HAARP  heater,  as  well  as  the  receiver  in 
Chistochina. 


boundary  of  the  Earth-ionosphere  waveguide,  it  is  not  clear  as 
to  whether  efficient  coupling  of  ELF  wave  energy  upward  (to 
DEMETER  altitudes)  is  indicative  of  efficient  coupling  to 
ground-based  sites  (and/or  vice  versa).  The  ELF  radiation 
directed  upward  from  the  heated  spot  into  the  magnetosphere 
may  be  significantly  different  from  that  directed  into  the 
Earth-ionosphere  waveguide.  Figure  4b  shows  a  schematic 
of  the  heated  region  and  illuminated  region,  with  respect  to 
the  ground  receiver  at  Chistochina.  The  vector  I  represents 
the  electrojet  current,  while  dl  is  a  differential  vector  length 
along  I.  Also  plotted  are  He  and  the  electric  and  magnetic 
field  vectors  observed  at  the  ground  receiver. 


5.  Calculation  of  the  ELF  Power  Radiated  by 
HAARP 


[ii]  To  estimate  the  total  ELF  power  radiated  as  a  result 
of  modulated  HE  heating  by  HAARP,  we  use  the  expression 
below  [Jasna  Ristic-DJurovic,  1993,  pp.  40-41]  to  calculate 
the  time  average  Poynting  flux  vector  (Sd),  using  local  cold 
plasma  parameters  (as  measured  on  DEMETER)  as  well  as 
the  value  of  the  total  wave  magnetic  field  Bq,  for  arbitrary 
values  of  wave  normal  angle  \|/: 


(^d) 


\l (tan\[/  -  pip2X)^+(l  +  pjX) 


COSVJ/) 


1  +  (1  +  tan^  v|/) 


(^2  sin2 

^  P{n^-S) 


12' 


(i) 


where  (...)  symbol  denotes  averaging  in  time  over  the 
duration  of  each  ELF/VLF  pulse,  where: 


,Ez  —  *5)^2^  sin \j/ cos \j/ 

' Ey  sin^  \j/ —  P)  ' Ey 


andX  = 


P 

[P  —  tP  sin^  \j/) 


D 

(2) 


The  refractive  index  n  as  well  as  the  values  of  the  auxiliary 
variables  D,  P,  and  S  are  obtained  from  Stix  [1962,  p.  10]. 
These  values  depend  on  Hg,  Hi,  the  electron  and  ion  plasma 
frequency  and  Gg,  Di,  the  electron  and  ion  gyrofrequency  as 
well  as  the  wave  angular  frequency  u;.  The  gyrofrequency 
values  are  calculated  using  a  DGRF/IGRF  Geomagnetic 
Field  model  with  Bq  =  40.804  piT  at  the  HAARP  geographic 
coordinates.  The  plasma  frequency  values  are  determined 
from  measurements  onboard  DEMETER  of  electron  and  ion 
(mostly  O^)  density  using  the  Langmuir  Probe  Instrument 
(ISL)  and  the  Plasma  Analyzer  Instrument  (lAP),  for 
an  electron  density  of  22^  ^  4750  cm“^  and  an  ion  density 
of  fii  ^  5352  cm“^.  The  refractive  index  n  also  depends  on 
the  wave  normal  angle. 

[12]  As  shown  in  Figure  2,  the  HAARP-generated  ELF 
signal  is  observed  on  DEMETER  at  a  distance  for  which  the 
subsatellite  point  on  the  ground  is  430  km  away  from 
HAARP.  The  signal  is  observed  during  a  30  second  period, 
after  which  the  DEMETER  receiver  is  shut  off  Therefore 
our  observations  are  representative  of  only  the  outer  part  of 
the  illuminated  region.  In  order  to  integrate  Sd  over  the  total 
area  and  obtain  an  estimate  for  the  radiated  ELF  power,  we 
must  assume  a  model  of  the  distribution  of  ELF  wave  power 
over  the  illuminated  area  B.  We  assume  that  the  variation  of 
Sd  is  symmetrical  around  the  magnetic  field  line  intersect¬ 
ing  HAARP  and  we  define  a  function  g(r)  which  represents 
the  radial  variation  of  the  Poynting  vector  with  r,  the  radial 
perpendicular  distance  from  the  HAARP  magnetic  field  line. 

We  define  a  function  Sm{r)  =  (r/2m)^+  (85km)^  as 

the  total  ray  path  from  the  heated  region  to  the  injection  point 
in  the  ionosphere  at  a  radial  distance  r,  for  a  ray  reflecting  m 
times  in  the  Earth-ionosphere  waveguide,  highlighted  in 
green  in  Figure  2  for  w  =  1 .  Assuming  that  the  power  varies 
as  were  ^  is  the  distance  to  the  source,  the  radiated  power 
Pt>  is: 


g(r)r  dr  d^ 


=  {^d) 


r\  n 

Ps  /  / 

^  m=l  i 


I'Krotil  -  CT,)" 


V  s„{r)  ) 


2 

dr 


(3) 


where  cp  is  the  angle  around  the  HAARP  magnetic  field  line, 
the  radial  distance  from  the  r^  is  the  radius  of  the  illumi¬ 
nated  region  C,  (ri  ^21  km).  It  should  be  noted  that  during 
the  30  seconds  of  observation  of  the  HAARP  signal, 
DEMETER  moves  a  distance  of  approximately  6.7  km/sec 
30  sec  =  200  km  over  the  illuminated  region.  Therefore  it 


4  of  6 


L16101 


PLATING  ET  AL.:  THE  HAARP-DEMETER  EXPERIMENT 


L16101 


Table  la.  Results  From  the  Caleulations  of  the  Estimated  ELF 
Power  Radiated  From  the  Eleetrojet  Virtual  Antenna  Into  the 
Magnetosphere,  Pd,  Using  Equations  (l)-(3)  Using  \|/  =  0°,  at 
Eaeh  of  the  Transmitted  Step  Frequeneies,  for  the  Distribution 
Model  Adopted 


/[Hzl 

SrY  Distribution 

(So),  pWW 

(Pd),  W 

118 

n/a 

n/a 

254 

9.4 

2.7 

367 

1.13 

0.32 

707 

17.5 

4 

798 

18.4 

3.4 

should  be  noted  that  the  observations  at  different  frequen¬ 
cies  occur  a  different  locations  inside  the  illuminated  region. 
We  assume  that  the  size  of  this  region  is  the  same  for  all 
frequencies,  but  the  change  in  location  of  the  spacecraft 
must  be  taken  into  account  when  calculating  g(r).  This  is 
achieved  introducing  the  parameter  r2(/),  the  distance  of 
the  satellite  footprint  to  HAARP  at  the  beginning  of  the  ELF 
pulse  of  frequency  /  (r2(/)  =  207  km  @118  Hz,  205  km  @ 
254  Hz,  197  km  @  367  Hz,  159  km  @  707  Hz  and  125  km 
@  798  Hz).  The  first  term  of  the  sum  on  the  right  hand  side 
of  (3)  represents  the  wave  energy  radiated  directly  upward 
from  the  heated  spot.  The  second  term  represents  the  energy 
leaked  upward  from  the  Earth-ionosphere  waveguide.  The 
quantity  cr^  represents  the  fraction  of  wave  power  transmit¬ 
ted  upward  into  the  ionosphere  as  the  waves  in  the  Earth- 
ionosphere  waveguide  reflect  from  the  upper  boundary.  We 
assume  that  Ot  is  constant  for  each  ray.  The  ionosphere  wave 
power  transmission  coefficient  cr^  at  the  time  of  the  observa¬ 
tion  is  unknown,  since  the  ionospheric  density  profile  in  the 
D-region  at  that  time  is  not  known.  Nevertheless  we  can  use 
(3)  to  determine  a  lower  bound  to  the  amount  of  power  that 
HAARP  is  radiating  into  the  waveguide,  if  we  choose  a  value 

of  that  maximizes  the  term  ^  cr^(l  — 

m=l 

is  found  that  this  term  is  maximum  for  cr^  approaching  1. 

The  factor  is  the  spread  of  power  as  the  rays  propagate  in 
the  magnetosphere  from  B  to  A  and  is  equal  to  the  surface 
ratio  A/B,  that  is  =  (0.69)“ \ 

[13]  Shown  in  Figure  4a  are  the  measurements  of  one  of 
the  components  of  the  horizontal  magnetic  field  at  Chisto- 
china,  which  we  denote  as  Bq.  If  we  assume  that  the 
electrojet  virtual  antenna  can  be  described  to  first  order  as 
a  dipole  with  current  I  in  the  East- West  direction  and  a 
moment  I  dl  (Refer  to  Figure  4b),  we  can  use  the  measure¬ 
ments  of  Be  from  the  ground  station  to  estimate  the  time 
average  Poynting  flux  (Sen)  of  the  wave  propagating 
towards  the  ground  and  the  power  radiated  towards  the 
ground  Pen  as: 


(Sc//)=x(E^  xHe)  = 


ho 


2\4{\+F^)(\+rY 


-Gr  and 


PCH  = 


J  (^ch)  •  dS  =  j j  {ScH)rl  sin  (\)  d(\)dd 

0,c^ 


(4) 


where  0  is  the  angle  around  the  electrojet  antenna  dipole  and 
ro  =  92  km,  as  shown  in  Figure  4b.  Chistochina  is  in  the  far 
field  of  the  modulated  electrojet  currents  at  frequencies 
greater  that  ^600  Hz,  and  in  the  near  field  for  lower 
frequencies.  At  these  lower  frequencies  the  near  field  terms 
must  be  subtracted  from  the  measured  fields  at  Chistochina 
before  it  can  be  used  for  power  calculations.  This  subtrac¬ 
tion  was  carried  out  assuming  that  the  modulated  electrojet 
currents  represent  a  dipole  current,  resulting  in  the  term  (1  + 
F^)  in  (4),  where  F  =  c/(ijj  rg).  The  magnetic  field  observed 
on  the  ground  receiver  at  Chistochina  is  actually  a  factor  of 
(I  +  R)  greater  than  the  actual  value  emitted  by  the  electrojet 
antenna,  due  to  the  reflection  on  the  ground,  where  R  is 
the  reflection  coefficient  of  the  ground.  For  dry  soil  (cr  ^ 
10“^  S/m  and  ^  4)  we  estimated  the  value  of  R  as  R  = 
Re[(^CH  —  1)/(^CH  +  1)]^  where  ^ch  is  the  refractive  index 
of  the  ground,  ncu  =  (£r  +  icr/(eo  @)^^^. 

6.  Results 

[14]  Results  of  the  power  calculations  are  summarized  in 
Tables  la  and  lb.  In  calculating  the  power,  ionospheric 
losses  were  neglected  since  these  are  generally  <1  dB  for 
the  frequencies  of  interest  here  [Helliwell,  1967].  Table  la 
show  the  results  for  the  power  calculations  using  the 
distribution  of  Sj)  along  the  region  B  explained  in 
section  5  by  equation  (3).  The  values  of  Sj)  are  obtained 
from  equations  (1)  and  (2)  using  =  0°.  Observed  values  of 
magnetic  field  on  the  ground  and  are  obtained  using 
the  N/S  and  E/W  oriented  antennas  at  Chistochina 
(Table  lb).  The  values  of  the  magnetic  field  for  each 
component  are  averaged  over  the  pulse  duration.  ^Sch  and 
Pen  are  the  Poynting  vector  and  radiated  power  at  each 
frequency,  respectively,  estimated  using  equation  (4).  From 
the  simple  model  used  to  determine  the  total  ELF  power 
radiated  by  HAARP,  we  arrive  at  results  which  resemble 
those  of  similar  previous  experiments  [James  et  al,  1990], 
with  a  minimum  estimated  value  of  radiated  power  of  5.9  W 
at  707  Hz.  Results  reported  in  James  et  al  [1990]  were 
performed  with  an  ionospheric  heater  located  at  a  similar 
latitude  as  HAARP,  the  MPAe  facility,  near  Tromso, 
Norway  (^69.6°  invariant  latitude),  radiating  with  an 
effective  HE  power  of  270  MW,  using  an  X  and  O  mode, 
AM  modulated  ELF/VLF  signal  that  consisted  of  a  series  of 
pulses  ranging  in  frequency  between  525  Hz  and  5925  Hz. 
The  DE  1  spacecraft,  which  observed  the  signals  during  this 
experiment,  was  located  at  an  altitude  of  ^11000  km,  ten 
times  the  altitude  of  DEMETER  in  the  present  experiment. 
Nevertheless  the  power  estimate  of  ^4W  is  comparable  to 
the  estimates  presented  in  this  report.  The  values  of  power 


Table  lb.  Results  of  Radiated  Power  Into  the  Earth-Ionosphere 
Waveguide,  Pen,  Obtained  From  Equation  (4)^ 


A  Hz 

W,  PT 

pT 

(ScH>,  pW/m" 

(Pch),  W 

118 

n/a 

n/a 

n/a 

n/a 

254 

1.54 

1.23 

50.7 

2.71 

367 

1.25 

0.87 

52.7 

2.82 

707 

0.92 

0.91 

78.8 

4.22 

798 

0.87 

0.78 

69.7 

3.73 

^Also 

included  are 

the  averaged  values  of  the  two 

magnetic  field 

components  measured  by  the  E/W  and  N/S  antennas  at  Chistochina. 
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estimated  using  DEMETER  wave  amplitudes  are  generally 
lower  than  the  estimated  power  based  on  ground  observa¬ 
tions.  In  particular,  the  power  estimated  using  DEMETER  at 
367  Hz  and  647  Hz,  reveals  a  strong  attenuation  that  is  not 
observed  on  the  ground  receiver.  The  pulses  at  frequencies 
below  the  local  proton  gyrofrequency  are  subject  to  mode 
conversion  from  the  right  handed  to  the  left  handed  mode, 
which  may  be  affecting  the  pulses  at  647  Hz  and  367  Hz.  It 
can  be  seen  that  there  is  a  modest  frequency  dependence  of 
the  radiated  power,  which  can  be  interpreted  as  an  increase 
with  frequency  of  the  efficiency  of  the  generation  of  ELF 
waves  by  modulation  of  the  electrojet  (in  space:  James  et  al 
[1984],  on  the  ground:  Rietveld  et  al  [1989]).  However  this 
dependence  is  less  than  the  f  dependence  which  would  be 
expected  from  a  dipole  antenna  of  fixed  length  in  free  space. 
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